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Abstract 

*^3 ' When time reversal is broken the viscosity tensor can have a non vanishing 

^L ■ odd part. In two dimensions, and only then, such odd viscosity is compatible 

j— with isotropy. Elementary and basic features of odd viscosity are examined by 

-2 \ considering solutions of the wave and Navier-Stokes equations for hypothetical 

O ■ fluids where the stress is dominated by odd viscosity. 

• i-H ■ 
C/3 , 
>v 

1 Introduction and Overview 

pg ' Normally, one associates viscosity with dissipation. However, as the viscosity is, 

in general, a tensor, this need not be the case since the antisymmetric part of a 

£j~ tensor is not associated with dissipation. We call the antisymmetric part odd. It 

must vanish, by Onsager relation, if time reversal holds. It must also vanish in 
three dimensions if the tensor is isotropic. But, in two dimensions odd viscosity is 

^ ' compatible with isotropy. 

It is conceivable that that odd viscosity does not vanish for many system where 
time reversal is broken either spontaneously or by external fields. But, I know of 
only two systems for which there are theoretical studies of the odd viscosity and 
>>> none for which it has been studied experimentally. In superfluid He 3 , time reversal 
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and isotropy can spontaneously break and the odd viscosity has three independent 

components ||. As far as I know there is no estimate for their magnitudes. In the 

two dimensional quantum Hall fluid time reversal is broken by an external magnetic 

field. In the case of a full Landau level the dissipative viscosity vanishes. The odd 

viscosity, rj a , has been calculated for non interacting electrons in M for the lowest 

Landau level. Using results of ||, the odd viscosity for any integral filling factor n 

is: 

n 2 pB 

V a (n) = Y e —, n = l,.... (1) 

e is the charge of the electron, c the velocity of light. It is an amusing coincidence 
of the cgs units that the kinematic viscosity of the electron gas, at integer fillings, 

8m,' (2) 




is close to n in cgs units. 

Significant odd viscosity can be responsible for odd properties and we shall il- 
lustrate these by considering the wave and Navier-Stokes equations for hypothetical 
media where the odd viscosity dominates. A good example for such a peculiar prop- 
erty is the following: Consider a small, slowly rotating circular cylinder (so that 
the Reynolds number is small) in an fluid which has both dissipating and and odd 
viscosity. The dissipative viscosity applies a torque which resists the rotation. This 
agrees with common intuition. The odd viscosity leads to radial pressure on the 
cylinder which is proportional to the rate of rotation (and the coefficient of odd 
viscosity, of course). Reversing the orientation of rotation, reverses the sign of this 
pressure. This response to the rotation is not particularly intuitive. 

Materials with significant odd viscosity, be they solids or liquids, can support 
non-dissipating chiral viscosity waves with a quadratic dispersion. The reflection 
from a boundary of these waves obeys a different rule than the reflection of acoustic 
waves. In addition, such materials can function as circular polarizers for ordinary 
acoustic waves. At the same time, viscosity waves can be quite elusive. In particular 
we shall show that there are no viscosity waves in isotropic and incompressible media. 
In particular, the Hall system, being isotropic and incompressible, does not support 
viscosity waves. 

A dimensioless number that gives a measure of the importance of viscosity rela- 
tive to elasticity is 

£=— (3) 

where v is the kinematic viscosity, uj the frequency, and c s the velocity of sound. 
Odd viscosity is always unimportant relative to elasticity at low frequencies. For the 
dissipating viscosity, and also for the odd viscosity in the Hall effect, the kinematic 
viscosity is of order 1 in cgs units. Then, since c s is typically of order 10 5 in cgs, e 
is of order unity for uj of order GHz. 

As we shall see, the generalized Navier-Stokes equation that allows for odd vis- 
cosity preserves the basic properties of fluid dynamics of the ordinary Navier-Stokes 
equation: Kelvin theorem and Bernoulli law generalize to non zero odd viscosity. 

2 Odd Viscosity 

Consider a hypothetical, homogeneous and ideal Newtonian fluid. The stress due to 
viscosity is: 

a ij = VijklUkl- (4) 

where Uij = | (wjj + Uj ti ) is the strain, iij k is the strain rate, r\ a (constant) viscosity 
tensor and a the stress. By general principles, |l|], the viscosity tensor, rjijj-i, is 
symmetric under i -«->• j, and k <r-> I. One can always write rj = r/ s + r\ where 
r] A is anti symmetric under {ij} *-* {kl} and r] S is symmetric under {ij} <-► {hi}. 
By Onsager relation || the anti-symmetric part is odd under time reversal and the 



symmetric part is even, e.g., with B an external magnetic field; 

V S (B)= V S (-B), V A (B) = - V A (-B). (5) 

So, for i] A ^ time reversal must be broken. One can ask if odd viscosity is not 
just a way of mascarading, say, an external magnetic field. The answer to this is 
no. This can be seen by a counting argument. A (constant) magnetic field has three 
components while the space of anti-symmetric tensors r\ A is 15 dimensional in three 
dimensions. 

2.1 Viscosity in Two Dimensions 

In two dimensions there is a natural basis for representing real 4-th rank tensor 
which are symmetric in pairsQ: 

»7 = S>«6 *° ® <A a, be {0,1, 3}. (6) 

Here a a are the Pauli matrices with <t° the unit matrix, a 2 , which is conventionally 
imaginary, is not used. In components: 

vw = Y.W}Ai ^J'e {0,1,3}. (7) 

The space of such tensors is nine dimensional, and its splits to a six dimensional 
even part and three dimensional odd part. The odd part takes the form 

/ = I>j,(cT a <g>C7 b -C7 a ®CT 6 ) Z,J6{0,1,3}. (8) 

To see what happens in isotropic media, recall that ia 2 is the generator of rotations 
in two dimension, and it anticommutes with a a for a = 1, 3 and commutes for a = 0. 
Hence a 2 <8> cr 2 commutes with a a <S> a b if a, b = 1,3 and a = b = 0. The isotropic 
symmetric part is two dimensional and is characterized by two viscosity coefficients: 

v s = v s (a 1 <g> a 3 + a 3 <8> a 1 ) + ( a ® a . (9) 

The odd isotropic part is one dimensional and is of the form 

v A = v a {a 1 ® a 3 -a 3 ® a 1 ). (10) 

The non-zero components are determined by 

r?H22 = 0, r] A u = - V A 22 = r] a . (11) 

It is sufficient that a two dimensional medium is invariant under rotation by 7r/4 for 
rj to have this form. 



lr This basis was suggested by the referee of this paper. 



As a consequence of this, isotropic two dimensional media are characterized, in 
general, by three coefficients of viscosity: two for the even part and the third for 
the odd viscosity if . This is in contrast with a claim made in || (p. 45-46), and in 
[|] §15, that isotropy alone implies that the viscosity tensor has two coefficients of 
viscosity. 

The stress a v associated with the viscosity of an isotropic medium, is 

<T% = 1] S (Uij + Uji) + C^ijUkk 

-2r] a (5 il 5 jl - Si2S j2 )ui2 + V a ($ii$j2 + 8 a 6ji)(un - itm)- (12) 

In an incompressible fluid iikk = and the stress becomes independent of (, the 
second viscosity coefficient. Because of this ( plays no role in the Navier Stokes 
equation for incompressible fluids. Incompressible and isotropic fluids in two di- 
mensions with broken time reversal are characterized by two viscosity coefficients 
one for the odd part, r] a , and one for the even part r] s . 

2.2 Viscosity in Three Dimensions 

In three dimensions isotropy implies that r\ A = 0. Three dimensional isotropic 
media are characterized by only two coefficients of viscosity, and both of these are 
associated with the even part of the viscosity tensor. Isotropic and incompresible 
media are characterized by a single dissipative viscosity coefficient. 

The vanishing of the odd viscosity for isotropic tensors can be seen by the follow- 
ing elementary argument, which I owe to L. Sadun. Vectors in three dimensions are 
associated with the 1 representation of the rotation group. 2-tensors are identified 
with 

1® 1 = 2©1©0 

representation of the rotation group. The vector representation 1 on the right hand 
side, is associated with pseudo vectors and so with the antisymmetric 2-tensors. 
The 2 © is the six dimensional representation associated with the symmetric 2- 
tensors. This shows that there is one isotropic 2-tensor in three dimensions, namely 
the identity. Continuing in this vein, tensors tijki with the symmetry % «-► j, and 
k <-> I are identified with 

(2©0) 2 = 4©3©3-2©l©2-0. (13) 

This shows that the space of isotropic 4-tensors with the symmetry % «-> j, and k ^ I 
is 2-dimensional. The isotropic 4-tensors t s make a two dimensional family, given 
by 

tijki — hSijSki + 1 2 (SikSji + SuSjk) , (14) 

this leaves nothing for the isotropic antisymmetric 4-tensor and so r\ A = 0. 

The number of independent components of the odd and even part of the viscosity 
tensor is given in the table. 



Dimension 


Even 


Odd 


general 


isotropic 


general 


isotropic 


2 


9 


2 


3 


1 


3 


21 


2 


15 






3 Viscosity Waves 

Consider a hypothetical ideal three dimensional medium where the stress tensor is 
dominated by the odd viscosity. For small oscillations the equation of motion is M : 



pui 



Wr 



(15) 



where p is the (mass) density, Uj the (cartesian) j-th component of the displacement 
(assumed small), and a v , the stress tensor, depends linearly on the strain rate. 
Linear elasticity is mathematically trivial in the sense that solving the PDE reduces 
to a problem in linear algebra. In a homogeneous systems p is a constant so one 
can, without loss, set p = 1 (by absorbing p in a.) For a medium with odd viscosity 
Newton equation gives a linear PDE which is second order in space and time: 



U: 



A 
Vijkl U l,jk- 



(16) 



This equation is not just an alternative way of describing the Lorentz force. The 
corresponding Newton-Lorentz equation: 



ii + B Au = 0, 



(17) 



is not a PDE. 

Newton equation for viscosity waves admits an integral and by choosing the 
constant of integration appropriately, reads: 



Ui = VijkiUijk- (18) 

This equation is first order in time and second order in space. Because r] A is anti- 
symmetric it has the same character as Schrodinger equation^. But it is a classical 
equation in the sense that the wave function is directly observable. 

Consider plane waves solutions to this equation with wave vector k. Let fi(k) 
be the (pseudo) vector with components that are quadratic forms in the coordinates 
kf 

Qi(k) = t ijk rjf apk k a k p . (19) 

Newton law for a plane wave propagating in the k direction can be written as: 



u = n(k) a u. 



(20) 



2 In real coordinates, with ip\ = Reijj and i\)i = Imip, Schrodinger equation for a free particle 
takes the form ipi = tijipj t kk- 



From this equation it is clear that the viscosity wave in three dimensions is chiral and 
circularly polarized in the plane perpendicular to the vector f2(k) and has quadratic 
dispersion. Recall that ordinary sounds waves have linear dispersion, but, rods and 
plates also allow for elastic modes with quadratic dispersion [|TJ. The equation for 
viscosity waves equation of motion, Eq. (|20|), is formally identical to the Landau 
Lifshitz equation for magnons fj 

3.1 Incompressible Media 

Consider an incompressible medium where V ■ u = 0. For a plane wave this means 
that ft is parallel to k, that is: 

nAk=0,^0 = ^M (21) 

For a viscosity wave propagating in the 1 direction 



n(k) = k 
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Vim | • (22) 

V<112. 



From this it follows that in an incompressible medium either r]^ u = for all i or 
u = identically. This leads to: 

Proposition 3.1 There are no odd viscosity waves in isotropic and incompressible 
fluids. In particular, there are no viscosity waves in a two dimensional quantum 
Hall fluid. 

This is so because in two dimensions isotropic fluids the non zero component of the 
isotropic odd viscosity are rf = —77^22 = ^2111 • For an isotropic fluid we can choose 
the 1 axis to coincide with the wave vector, and then incompressibility says that 
either r] a = r\^ XY y = or the wave has zero amplitude. [] In three dimensions, there 
are no viscosity waves because the odd viscosity vanishes. 

3.2 The Energy Flux of Viscosity Waves 

In the case rf = one can define a conserved current associated with conservation 
of energy. The kinetic energy density: 

£ = \{Reu)\ (23) 



satisfies a conservation law: 



d t S = diJi, (24) 



3 I thank Dr. M. Milovanovic for reminding me of this fact. 

4 The result, and proof, works also for a two dimensional planar medium embedded in three 



dimensions. 



where the energy flux J is 

Ji = vtjki ^ (uj) Re{ii kl ) = a^Reiiij). (25) 

It follows that odd viscosity carries energy flux. There is a generalization of this 
energy flux to media where the stress also has elastic part. 

3.3 Dissipative Shear Waves 

It is instructive to contrast viscosity waves with the strongly dissipative shear waves 
[H]. For shear waves rj s ^ and r\ A = 0, and the dispersion curves are pure imaginary 

ito = -r] S k 2 , rj s ,uj>0 lmk>0. (26) 

Shear waves decay rapidly, by a factor of about exp(27r) rs 540 in one period Q. 
If one now consider a mixed situation where 77 and r/ A are of the same magnitude 
then the attenuation is much weaker — of order exp(27rtan(7r/8)) ~ 13. 

4 Scattering of Viscosity Waves 

Using the rules of geometric acoustics it is straightforward to find the laws of reflec- 
tion and refraction of an acoustic wave (in a medium where rj = 0) from an interface 
with an odd viscous medium (where the stress is dominated by r/ A ). If we let y 
be the axis separating the two media, then fc 2 is conserved and so is the frequency 
to. Let 6 be the angle of incidence, cos(6>) = x ■ k\, of an acoustic wave with sound 
velocity c s and the angle of the transmitted viscosity wave, cos(0) = x ■ k t . Then 
a little geometry shows that 

sm 2 (f) = esm 2 e, e = ^. (27) 

pel 

We see that e^ 1 ! 2 plays the role of the index of reflection. Since e is small at low 
frequencies, the index of reflection is large, and the transmitted beam is essentially 
normal to the interface. Since the index depends on u the reflection is dispersive. 

4.1 Reflection from Empty Space 

Consider the reflection of ideal viscosity waves in two dimensional isotropic and 
compressible medium from an interface with a vacuum which we take to be the line 
x = 0. A viscosity waves is 

u(z) = f . 1 z, z = expi(k ■ x — u>t), (28) 

where u = r] a k 2 . The waves is circularly polarized in the plane and changes polar- 
ization from longitudinal to transversal. The boundary conditions at x = are 

o-n = a 12 = 0. (29) 



For this wave: 

ctii(u) = —771112 k, cr 12 (u) = — M71112 k; k = k 1 + ik 2 , k = ki — ik 2 - (30) 

Let k = (ki,k2) be the incident wave, k = (— ki,k 2 ) be the reflected wave; and 
z = exp i (k • x — uot) . u(z) is the incident wave (from the left) and u(z) the reflected 
wave and R(k) the reflection amplitude. The wave on the left is: 

u(z) + R{k)vL(z), x<0. (31) 

Substituting in the boundary conditions, using Eq. (|5D|) and linearity, one finds 

R(k) = t. (32) 

The phase of the reflected wave has an interesting dependence on the direction of 
incidence: It goes from 1 for normal reflection to —1 for grazing reflection. This 
is quite unlike the reflection of say longitudinal sound waves in an incompressible 
liquid where the reflection amplitude is independent of the direction of incidence. 

4.2 Odd Viscosity as Polarizer 

Consider a wave propagating in the x direction, of a three dimensional incompressible 
medium, so that to the left there is an isotropic elastic medium with velocities of 
sound c t and q. To the right is a viscous (possibly incompressible) medium. The 
y — z plane separates the two media. Incompressibility says that f2(k) points in the 
x direction. Let 

e + = (0,1, z), e_ = (0,l,-z), (33) 

denote the two basic vectors of circular polarization. Consider scattering of a trans- 
verse wave with positive chirality from the y—z interface. The incoming and reflected 
waves in the elastic medium are 

e + (expik(x — c t t)+rexpik(—x — c t t)), x < 0. (34) 

The transmitted wave (necessarily with positive chirality) is: 

t e + exp i(kx — out) , x > 0. (35) 

Where uj = kc t = |f2(k)|. The basic boundary conditions matches ah on the left 
with afj on the right. One checks that af x = a^ = and that ia^ 2 = a^ 3 . ia^ 2 = o'is 
is automatically satisfied ( confirming the anzatz that there is no reflected wave with 
flipped chirality). The remaining equation, crf 2 = af 2 i gives 

1 — r = et, e= ^— > 0. (36) 



Continuity of the wave u on the boundary gives two equations but only one new: 
1 + r = t. From these I finally get: 

, 2 1-e T] 2113 U 

1 + e 1 + e pet 

Note that \r\ < 1 as it should: The reflected wave always has smaller amplitude 
than the incident wave. However, for t one gets an unusual behavior: < t < 2 and 
the amplitude t is maximal for e = when the amplitude of the transmitted wave 
is twice as large as that of the incident wave! This surprise is mitigated when one 
notices that for the energy flux in the x direction, Ji, the dependence on e is (1 f £)2 , 
which gives maximal flux at e = 1 as one expects, since in this case r = 0. The wave 
is perfectly transmitted. 

Consider now the scattering of a viscosity wave so that the incident wave has 
negative chirality. The most general transverse wave in the elastic medium with the 
given incident wave is: 

re + exp ik(—x — q£) + e_ ( exp ik(x — q£) + f exp ik(—x — c t t) J , (38) 

while the wave in the viscous medium is as before. Matching erf- with a^ one finds 
that as before a^ = crf^ = and ia^ 2 = cfl- Writing out icrf 2 = crfs gives: f = 1. 
The remaining equation, <rf 2 = a^ 2 , gives 

r = et. (39) 

This is all that follows from the continuity of the three stress components. Imposing 
continuity of the wave on the boundary gives two new equations: 

r = t, l+f = 0. (40) 

The second equation is in conflict with the equation for the continuity of the stress^. 
The way out is not to require continuity of the wave, but instead continuity of the 
energy flux. This holds if the flux vanishes at the surface and sets r = t = 0. 
The incident wave with the wrong polarization is totally reflected. We see that 
anisotropic incompressible media with odd viscosity act as circular polarizers. 

5 Navier Stokes Equation 

Consider the Navier Stokes equation, for homogeneous (p = 1), isotropic and in- 
compressible fluid in two dimensions-the one dimension where odd viscosity is com- 
patible with isotropy. 



5 This is the generic situation in scattering of acoustic and viscosity waves. That is, in general 
imposing continuity of the wave across the interface gives an overdetermined system. In this 
sense, the previous example is exceptional in that both the stress and the wave turned out to be 
continuous. 



With broken time reversal, the general Navier-Stokes equation for the pressure p 
and the velocity field v, is the obvious generalization of the standard Navier-Stokes 
equation ||: 

<9 t v+ (v V)v = -Vp + n s Av-n a Av*; V • v = 0. (41) 

The dual is defined, as usual, by 

v* = e ijVj . (42) 

Incompressibility is implemented by introducing a stream function and the equations 
can be alternatively written as four equations for the four fields (p, ip, v): 

dtv + (v • V)v = - V(p - V a + »] s Av; v=(V^)*, ^ = -AV>. (43) 
Taking curl the equation for the vorticity is 

dtt + (v • V) e = V s A£. (44) 

Vorticity is be generated only by the symmetric (dissipative) part of the viscosity 
tensor. The odd viscosity does not generate vorticity. 

5.1 Bernoulli Law 

Consider incompressible fluid in two dimensions, with n s = 0, in steady state. The 
equation of motion is: 

( v .V)v = -V(p-77°0. (45) 

Since 

- V(v 2 ) = (v • V)v + v x (V x v) = -V(p - r] a + v x (V x v), (46) 

Li 

one gets, by integrating along a streamline, that 

- v 2 + p - 77°£ = const. (47) 

This looks at first like an interesting generalization of Bernoulli's law to odd viscosity 
with the pleasant feature that the vorticity (weighted by the odd viscosity) plays a 
role of (signed) kinetic energy. However, as vorticity is conserved along a streamline, 
it is actually precisely equivalent to Bernoulli's law. 

5.2 Stokes Equation 

The limit of small Reynolds number is described by Stokes equation. This is the case 
where normally viscosity dominates. Stokes equation is linear, so this is also the easy 
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limit. In the present context, with broken time reversal, and for two dimensional 
isotropic and incompressible fluid: 

d t v = -V(p-ri a £)+ri s Av; v = (WO*. (48) 

Taking curl and div we get: 

c^ = r/ s A£, = -A(p-77 a £). (49) 

In a steady state, if rj ^ 0, both the pressure and the vorticity must be harmonic 
functions and if> bi-harmonic. 

5.3 A Rotating Disc 

Consider a rigid disc which is rotating with constant (and small) angular velocity 
Q in a fluid subject to no-slip boundary conditions. Stokes equation can be solved 
explicitly in this case [^ ||. I will show that while the ordinary viscosity applies a 
torque that resists the rotation of the disc, the odd viscosity leads to a pressure on 
the disc which is proportional to the angular velocity. The pressure can be either 
positive or negative, depending on the sense of rotation. 
Written in complex notation Stokes equation read: 

8(p + ir]€) = 0, r] = r] s + ir] a . (50) 

Since if) is bi-harmonic and real 



— 2iip = za{z) — za(z) + b(z) — b(z). (51) 

Let v = Vi + iv 2 . Using 

v = -2idip, f = -2idv. (52) 

one finds 

v = -a(z) + za'(z) -¥{z), £ = 2i(a'(z) -a'(z)). (53) 

Since the velocities at infinity are finite, one must have[] 

00 a ■ °° b- 

a(z) = Y: % b(z) = b\ogz+ £ -j. (54) 

j=o z i=-i z 

On the surface of the circle z z = 1 and one can trade z for an inverse power of z, 
so: 

itlz = - f (j ~ 1)& f 1 + Qj -£(;- 2)5^ - h, \z\ = l. (55) 
3=0 z i=3 

Equating coefficients we get 

iQ = -b. (56) 

6 One can actually allow <z_i real. One can not allow a log term in a(z) if one is interested in 
continuous solutions for the velocity field. 

11 



All other coefficients are zero. That is 

v = ~ij, ^ = 0. (57) 

It then also follows that p is a constant. Using the velocity field one finds the strain 
on the surface of the circle: 

wn = -ftsin(20), U12 = ttcos(29), u 22 = ttsm(29). (58) 

Here x — \z\ cos 6. Let us write the stress as a = a s + o a where the first piece is due 
to the dissipative viscosity and the second is due to the odd viscosity. Then 

o£ = 2i7'uy. (59) 

And 

a a xl = 2r] a n cos(2fl), a a l2 = 2r, a Vt sin(20), a a lx = -2rj a Q cos(2fl). (60) 

Let n and t denote the unit vectors normal and tangent to the disc. The a calcula- 
tions gives 

a s nn = 2 V a n, al = 2 V s n. (61) 

This says that the symmetric viscosity resists the rotation by applying a torque 
on the rotating circle. The odd viscosity applies no torque, but instead, a normal 
pressure on the circle, which is proportional to the rate of rotation. 

Remark: One could have also asked what is the effect of odd viscosity on the drag 
and lift (the Magnus force). Unfortunately, for an incompressible, viscous fluid in two 
dimensions the question is moot: Stokes equation in two dimension does not admit 
steady state solutions that describe a moving disc with no slip boundary conditions 
in a fluid that is at rest at infinity. Two dimensions behave like one dimension and 
unlike three dimensions. This feature of the Stokes equation is known in classical 
fluid mechanics as the Whitehead paradox. 
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